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ABSTRACT 

We present a set of new characteristic frequencies associated with accretion disks 
around compact objects. These frequencies arise from persistent rotating patterns in 
the disk that are finite in radial extent and driven purely by the gravity of the central 
body. Their existence depends on general relativistic corrections to orbital motion and, 
if observed, could be used to probe the strong gravity region around a black hole. We 
also discuss a possible connection to the puzzle of quasi-periodic oscillations. 
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1 INTRODUCTION 

Timing observations of accreting X-ray binary systems have 
revealed luminosity modulation at a number of charac- 
teristic frequencies. Phenomenology of these quasi-periodic 
oscillations (QPOs) is quite rich. For a detailed review, 
see [McClintock & Rcmillard (200!^ and references therein. 
Some of the features are rather puzzling, such as stability 
of high frequency QPOs in black hole binaries and that in 
some systems they appear in pairs at 3:2 frequency ratio. 

QPOs in black hole systems are thought to arise from 
physical processes in accretion disks. Depending on where 
the oscillations reside, one can roughly divide models for 
QPOs in accretion disks into two classes: local and global. 
Local models tie down the oscillation frequency to a particu- 
lar place in the disk (like an edge or a hot spot). In this case, 
the question of what determines that place has to be an- 
swered. One line of argument is th at the location of the hot 
spot is de termined by a resonance ramowic^^^Ouzniay 
^001; Tor ok et akE oOS: Kluzniak 2005i: lAbramowicz^l200.'^^ 
This model has an attractive feature that the observed 3:2 
frequency ratio can be explained by non-linear mode locking. 
A hot spot can give rise to luminosity variation, for example , 
due to Dopplcr beaming dSchnittman fc Bertschingeill20o3 : 
[Schnittman 2005). However, a potential difficulty is to have 
a hot spot which is sufficiently bright. Achieving sufficient 
luminosity variation seems less problematic in global mod- 
els, in which modes occupy a larger region of the disk. Lin- 
ear perturbation a nalysis of t he accretion disk in diskosois- 
mology approach dWae oner et al. 2001; Perez ct al . 1997; 
ISilberdeit et alJl200ll : lOrtega-Rodn'guez et al.ll2002l) natu- 
rally solves the issue of spatial and frequency localization 
of modes. The 3:2 frequency ratio would be accidental for 
two fundamental diskoseismic modes, but it could arise from 
higher azimuthal gr-modes which are nearly harmonic. 
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In this paper, we describe a set of new characteris- 
tic frequencies which might be present in accretion disks 
around compact objects. To the best of our knowledge, they 
have remained unnoticed in the literature. These frequen- 
cies arise from rotating patterns in the disk which are quasi- 
stationary, finite in radial extent, and driven purely by grav- 
ity of the central body. We neglect self-gravity and the hy- 
drodynamics of the accreting matter. The main idea is simi- 
lar to the notion of density waves that give rise to the spiral 
structure in galaxies ( Lindbl ad 1963: L in fc ShrLl964 . 1966). 
although these patterns depend on general relativity rather 
than a distributed matter source for their existence. 

While it is tempting to identify the frequencies of these 
patterns with the source of QPOs, we cannot claim to have 
a complete model. The issues of how they are excited, how 
they translate to X-ray luminosity variation, and effects of 
pressure and viscosity need to be investigated in more de- 
tail. We will return to these points with some plausibility 
arguments in Section 2] 



2 ACCRETION DISK KINEMATICS 

A test particle in a circular equatorial orbit around a Kerr 
black hole has an orbital frequency (.Bardeen et aL.1972i) 

1 



with respect to Boyer-Lindquist time t, where r is the orbit 
radius and a is the dimensionless black hole spin parameter 
(a = cJ/GM^). We work in dimensionless units scaled by 
the black hole mass M (i.e., distances measured in units of 
GM/(? , times measured in units of GM/(? , etc.), and will 
further set G = c = 1. Here and later we will assume that 
particles co-rotate with the black hole. 

If perturbed from the circular orbit, the particle will 
undergo radial and perpendicular o scillations with epicy clic 
frequencies k and Q,±_ respectively iOkazaki et al.l[l987^ 
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Figure 1. Orbital frequency Q, radial epicyclic frequency k, and 
precession frequency a; = f2 — 2k of a 2:1 orbit in an accretion disk 
around a Kerr black hole with M = 10 Mq and a = 1/2. Preces- 
sion frequency exhibits a shallow negative minimum at r» 9.64. 
The radial extent of the rotating pattern with 20% deviation in 
frequency is marked by vertical bars. 



and 



^ 6 ^ 8a 



r3/2 
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The factors multiplying Q, on the right-hand sides of these 
expressions are general relativistic corrections. They are ab- 
sent in Keplerian mechanics, where both epicyclic and or- 
bital frequencies are all the same (k = Q.± = Q). The radial 
dependences of orbital frequency Q. and radial epicyclic fre- 
quency « for a typical rotating black hole are illustrated in 
Figure Circular orbits close to a black hole are unstable; 
the innermost stable circular orbit (ISCO) is located where 
vanishes. 

If orbital and epicyclic frequencies are the same, as they 
are for a Keplerian potential, the orbits are closed. However, 
if the potential deviates from 1/r (either because of general 
relativity corrections, as in our case, or due to a distributed 
matter source, as happens in galaxies), the two frequencies 
will in general be different, and the orbits will precess. The 
condition for an orbit to close in a frame rotating with fre- 
quency u is for the orbital and epicyclic frequencies to be 
commensurate, m{Q, — uj) = nn, which gives the precession 
frequency 

Tl 

10 = 0, K. (4) 

m 

The integers n and m determine the shape of the precessing 
orbit, and from here on, we will use the abbreviation n:m 
to refer to their values. Figure |5| shows the shape of 1:2 and 
2:1 orbits, which are representative of the deformed and the 
self-intersecting orbit classes. 

In general, precession frequency ui depends on r, and 
any pattern initially present will shear away as the disk ro- 
tates. However, if lj is approximately constant over some 
portion of the disk, collective orbit precession can lead to 



a nearly rigid pattern rotation. One example of this is 
the spiral structure in gala xies caused by the 1:2 mode 
l|Lindbladlll963 : iLin fc Shul [T964. 1966). The 1:1 mode in 
Kepl erian disks gives rise to a static one-armed spiral pat- 
tern jKatdll983l), which is seen in numerical simulations as 
well jHavasaki fc Okazaki"2005') . An attempt has been made 
to tra p the 1:1 mode in the region of the strong gravity LKatol 
Il990l) . but trapping depe nds strongly on the p ressure distri- 
bution within the disk (IKato fc Honmalfl99 j) . 



3 PERSISTENT PATTERNS IN THE DISK 

Our key observation is that for n > m, the precession fre- 
quency ui{r) develops a very shallow minimum at a radius 
r = r*, as illustrated in Figure^for the 2:1 orbit. Collective 
excitation of particles on orbits precessing at the same rate 
would lead to a pattern occupying a sizable portion of the 
disk around r* and rotating with little shear at a frequency 
ujp = aj(r»). Somewhat unusual are the facts that the pat- 
tern is counter-rotating and that the orbit closes in several 
rotations rather than a single one. 

In a real accretion disk, the collective particle motion 
would have to be excited by some dynamical mechanism. 
It might be complicated and could require numerical simu- 
lations of the disk to fully understand the driving process. 
In the present paper, we will be content with studying the 
kinematics only. We will set off a collective mode by select- 
ing appropriate initial conditions (as described below) and 
follow the pattern evolution by tracing the motion of indi- 
vidual particles making up the disk. The purpose is to find 
out if there is a long-lived pattern that survives the differ- 
ential rotation. 

If perturbed from a circular orbit at r = ro by a small 
displacement e in the radial direction, the trajectory of a 
test particle (to first order in e) is 

r{t) = ro -I- e sin(Kt -I- x) , (5) 
20, 

<f){t) = if + Q.t ^ T£cos(fi:t 4-x), (6) 

nr 

where tp and x are initial orbital and epicyclic phases, and 
1 



T = 



3 I 2a 
r 7372 
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is a relativistic correction factor (which, however, changes 
little in the region of the disk we are interested in). We 
populate the disk by spreading A'^ particles uniformly on a 
n:m orbit, with initial phases of a fc* particle 

27rm 



2?™ 



Xk = 



N 



fc, 



(8) 



while stacking the orbits in a radial direction at an angle a 
by giving the orbit located at rj a phase offset 



(9) 



This particle distribution leads to a spiral structure in the 
disk. Figure |3] shows the surface density contrast (smoothed 
with a Gaussian kernel) for patterns obtained by distribut- 
ing particles on stacked 1:2 (left) and 2:1 (right) orbits. The 
three rows of Figure|3show a time-lapse sequence of pattern 
evolution. The top row shows the initial conditions, and the 
second and the third rows show patterns at t = T2:i and 
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1:2 



2:1 




tion of the 2:1 pattern is markedly different. Signs of shear 
are clearly seen after a single rotation. However, even af- 
ter twenty rotations, there is still a pattern present around 
r, ~ 9.6 (which is exactly where the minimum of uj2:i oc- 
curs). As this time span corresponds to almost 50 orbital 
rotations at r*, the pattern is remarkably persistent. 




Figure 2. Sample closed orbits. 1:2 orbit (left) corresponds to 
2 epicycles per 1 rotation around a central body, while 2:1 orbit 
(right) corresponds to 1 epicycle per 2 rotations. 
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Figure 3. Evolution of density patterns obtained by distributing 
particles on stacked 1:2 (left) and 2:1 (right) orbits. Top row shows 
initial configuration, middle row - after one period of rotation of 
2:1 pattern elapsed, bottom row - after twenty periods. Density 
contrast has been enhanced by histogram equalization. 



t — 20T2;i correspondingly. T2:i — 2tt/uj2:i denotes a period 
of rotation of the 2:1 pattern. 

The frequency of the 1:2 orbit precession depends mono- 
tonically on the radius, so one expects differential rotation to 
destroy the pattern. Indeed, at t = T2:i, the spiral is seen to 
wind up, and by t = 20r2:i, it is wound up so tightly that 
the smoothing removes all traces of structure. The evolu- 



4 DISCUSSION 

In the last section, we have shown that an accretion disk 
around a compact object can support persistent rotating 
patterns due to the collective excitations of particles in the 
disk. Their existence depends on general relativity effects 
and is sensitive to the parameters of the central body but 
not to the accretion rate. All the frequencies in the prob- 
lem scale inversely proportionally to the central body mass. 
In addition, persistent pattern frequencies depend on the 
spin parameter. Figure 2] shows the rotation frequencies of 
the three lowest-order persistent patterns (2:1, 3:1, and 3:2) 
for a lOAfQ-mass black hole as the spin is varied. This de- 
pendence in principle could be used to measure the mass 
and spin of the central object, provided that the frequencies 
of two distinct modes are observed and identified correctly. 
One should note, though, that for multi-armed patterns (for 
example the 3:2 pattern which has two arms) modulation 
frequency could be a multiple of the rotation frequency. 

Several different persistent patterns could coexist in the 
accretion disk; however, it is likely that the lowest-order ones 
are strongly selected based on geometrical considerations. 
Surface density modulation of the 2:1 pattern is second- 
order in particle displacement e, while the 3:1 and 3:2 pat- 
terns are third-order. The cancellation of lower-order terms 
is directly caused by the multiple-fold geometry of the or- 
bits with n > 1. Being higher-order could explain why these 
persistent patt erns are not apparent in the linear perturba- 
tion analysis oflWa goner et al ]f200l'). Despite second-order 
scaling, the 2:1 pattern in Figure |3| (lower right) shows 1% 
amplitude of the surface density modulation for moderate 
displacement values (e/r, « 0.022). The radial extent of the 
2:1 pattern (Ar ~ 4) also appears to be wider than that of a 
fundamental g-mode (the width of which is proportional to 
cV^ and is estimated as Ar ~ 1 bv lWagoner et all (1200111 ). 

In this paper, we focused on kinematics and neglected 
particle interactions and the hydrodynamics of the disk. The 
extent to which this approximation is justified should be 
further investigated. Of critical importance for the model is 
understanding the excitation mechanism. It is possible for 
the spiral waves to be driven from the outer edge of the disk 
JSpruitlll98^ : ISpruit et al.lll987ll : however, whether that is 
sufficient to cause persistent patterns to appear remains to 
be seen. Both questions could be answered by turning to 
numerical simulations of the accretion disk hydrodynamics. 
However, that is a much more complicated problem, and we 
feel that it is beyond the scope of this paper, the intent of 
which is merely to point out the existence of new character- 
istic frequencies in the disk. 

It is plausible that the characteristic frequency of the 
collective motion will manifest itself in X-ray luminosity 
variation, but the exact mechanism responsible for the mod- 
ulation is not clear to us. Density pattern in the accretion 
disk need not be a direct cause. Particles weaving in and out 
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Figure 4. Persistent pattern frequencies of three lowest-order 
radial modes (2:1, 3:1, and 3:2) for a lOMQ-mass black hole as a 
function of black hole spin parameter a. 
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Figure 5. Persistent pattern frequencies of three lowest-order 
transverse modes (2:1, 3:1, and 3:2) for a lOMQ-mass black hole 
as a function of black hole spin parameter a. 



on self-intersecting orbits could lead to efficient gas heat- 
ing, possibly due to shock formation, and create a temper- 
ature pattern in the accretion disk (in a sense, an exten- 
sive "hot spot"). The picture of temperature modulation 
of the disk causing X-ray luminosity variations is not en- 
tirely satisfactory as quasi-periodic oscillations are seen pri- 
marily in the hard non-thermal component of the emission 
jMcChntock fc Remillardll2003h . That could indicate that 
the quasi-periodic emissio n is coming from a coronal region 
rather than from a disk (iTitar chuk fc Fioritd ^2004^ . It is 
possible that the transfer mechanism might involve a mag- 
netic field threading the disk (Blandford, unpublished). 

So far we have been talking about patterns arising from 
radial oscillations. It is worth mentioning that a similar thing 
could happen for transverse oscillations as well. The preces- 
sion frequency uj±^ = 0,— ^ fix also has a minimum if n > m. 
However, the minimum lies inside an innermost stable circu- 



lar orbit unless the black hole is spinning rapidly (a > 0.8 for 
2:1 orbit). The frequencies of the three lowest-order trans- 
verse modes are shown in Figure |S] Transverse particle ex- 
citations would lead to a corrugated accretion disk rather 
than a surface density pattern. 

To summarize, we have found a set of new characteristic 
frequencies associated with accretion disks around compact 
objects. Although many questions remain, it might be in- 
teresting to pursue this idea further and see if it could lead 
to a model of quasi-periodic oscillations in X-ray binaries. 
In particular, the numerical values of our characteristic fre- 
quencies and their independence of the accretion rate sug- 
gest an application to high-frequency QPOs in black hole 
binaries. 
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